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ABSTRACT 

The Infrared Array Camera (IRAC) on the Spitzer Space Telescope has revealed that a number 
of high-mass protostars are associated with extended mid-infrared emission, particularly promi- 
nent at 4.5-/im. These are called "Green Fuzzy" emission or "Extended Green Objects". We 
present color analysis of this emission toward six nearby (d=2-3 kpc) well-studied high-mass 
protostars and three candidate high-mass protostars identified with the Spitzer GLIMPSE sur- 
vey. In our color-color diagrams most of the sources show a positive correlation between the 
[3.6]-[4.5] and [3.5]-[5.8] colors along the extinction vector in all or part of the region. We com- 
pare the colors with those of scattered continuum associated with the low-mass protostar L 1527, 
modeled scattered continuum in cavities, shocked emission associated with low-mass protostars, 
modeled H 2 emission for thermal and fluorescent cases, and modeled PAH emission. Of the 
emission mechanisms discussed above, scattered continuum provides the simplest explanation for 
the observed linear correlation. In this case, the color variation within each object is attributed 
to different foreground extinctions at different positions. Alternative possible emission mecha- 
nisms to explain this correlation may be a combination of thermal and fluorescent H2 emission in 
shocks, and a combination of scattered continuum and thermal H 2 emission, but detailed models 
or spectroscopic follow-up are required to further investigate this possibility. Our color-color di- 
agrams also show possible contributions from PAHs in two objects. However, none of our sample 
show clear evidence for PAH emission directly associated with the high-mass protostars, several 
of which should be associated with ionizing radiation. This suggests that those protostars are 
heavily embedded even at mid-infrared wavelengths. 

Subject headings: scattering — infrared: ISM 



1. Introduction 

The UV radiation and outflows associated with 
high-mass protostars hold keys for understand- 
ing high mass star formation and the formation 
of clusters. It has been debated for many years 
whether UV radiation or stellar wind could stop 
mass accretion, preventing the formation of high- 
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mass stars (see Beuther et al.|2007 for review) . Of 
the scenarios proposed to solve this issue, disk ac- 



cretion is the most promising (see Cesaroni et al 



2007 for review). Furthermore, high-mass pro- 



tostars are often associated with clusters of low- 
mass protostars. Energetic outflows and intense 
UV radiation from high-mass protostars (or high- 
mass young stars) could affect their formation 
and/or evolution, by destroying circumstellar ma- 
terial and/or triggering new generations of star 
formation (see Stahler and Paha 2005, for review). 

The Infrared Array Camera (IRAC) on the 
Spitzer Space Telescope has provided excellent ca- 
pabilities for studying outflows and the effects of 
UV fields associated with star forming regions. 



Observations have shown the presence of extended 
infrared emission toward a variety of protostars 
and star forming regions. This emission is often at- 
tributed to shocks associated with outflows, scat- 
tered continuum in the outflow cavities, or poly- 
cyclic aromatic hydrocarbons (PAH) excited by 
UV radiation. Such emission is expected to cover 
more than one or all filter bands (3.6, 4.5, 5.8, and 
8.0 /im) of IRAC (e.g. |Reach et aT 



et~aT1[2006l [Tobin et al.||2007[ |2008| |Neufe 



2009 



2006 




Smith 


Neufe 


c 


et al. 



marize the bright H2 lines and atomic/ionic lines 
which can be observed in IRAC bands. These are 
often shown with conventional three-color images 
with different filters (in many cases blue, green, 
and red for 3.6, 4.5 and 8.0 /im, respectively; e.g. 



Noriega-Crespo et al. |2004 Marston et al.| 2004 


Rathborne et al.||2005[ |Smith et al.||2006[ | 
et al. 2007; Kumar and Grave 2007; Shephcrc 
20071 |Tobin et al.l [20071 [20081 IQiu et all 


Araya 
= eTaT 
20081 


Teixeira et al.|2008| Cyganowski et al.|2008( 


2009 


2011 
Zhan 


[Chambers et al. 2009 Morales et al. 


2009 


g and Wang 2009[ Simpson et al.||2009 


Var- 


ricatt||2011). 



This three-color method provides useful diag- 
nostics for the nature of these sources. The PAH 
emission is brighter at 8.0-/xm than the other 



bands (e.g., Reach et al. 2006 ). Hence it appears in 



red with the above color combination (e.g., Shep- 



herd et al.|2007[|Kumar and Grave|2007l|Qiu et al 



2008). Other extended emission tends to appear 
"green", while stars are usually more "blue" as 
the flux is larger at shorter wavelengths. Studies 
of low-mass protostars show that emission from 
both shocks and scattered continuum appear green 
in three colors images (see, e.g., Noriega-Crespo 



et al. 


2004 


Smith et al.||2006 Tobin et al. 2007 , 


2008 




Teixeira et al. 


2008 


). Although the adjust- 



ment of color contrast has been rather arbitrary 
in individual publications, such identification of 
emission mechanisms seems to work fairly well (see 
also Section 2 for details). 

A number of high-mass star forming regions are 
known to be associated with such "green" emission 
in the three color image with blue, green, and red 
for 3.6, 4.5, and 8.0 /xm (e.g. |Marston et al.||2004 
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Rathborne et al.| [20051 jSmith et al.| [20061 l Ar aya| 
et al.||2007[|Shepherd et al.||2007[|Qiu et al.||2008[ 



Morales et all [20091 |Simpson et Tall [20091 |Varri 
catt|2011 ) . Candidates for high-mass star forming 



regions identified in this manner are called "Ex- 



tended Green Objects" (e.g., Cyganowski et al 



2008) or "Green Fuzzies" (e.g., [Chambers et al. 
2009[|De Buizer and Vacca||2010| ). H 2 and/or CO 
emission in outflow shocks are often regarded as 
the primary mechanisms to explain such emission 
(|Rathborne et al.||2005[ |Araya et al.||2007[ |Shep 



herd et al.||2007[ Qiu et al.||2008[ |Morales et aT 



2009; Cygan owski et a l. 2008, 2009, [20TTJ [Cham- 



bers et al. 2009). Observationally, it is in many 
cases based on comparisons with the studies of 
the well known massive protostellar outflow DR 21 



dSmith et al.|2006l |Davis et al.|2007[ ) and the low- 



mass protostellar outflow in the HH 46 / 47 system 



Davis et al. 


( 


2007 


), and 



(2004) show that the 4.5-/xm emission in these ob- 
jects have a morphology very similar to H2 emis- 
sion at 2.12 /xm, indicating that the "green" emis- 
sion in these objects are associated with molecular 
shocks. The same trend has also been observed 
for several features observed in high-mass star 
forming regions (Qiu et al. 2008[ ). Several Green 
Fuzzies are associated with molecular outflows ob- 



served at millimeter wavelengths (e.g., Shepherd 
et al I [2007] IQiu et al.| [20081 |Cyganowski et al. 
2011), Class I methanol masers, i.e., a tracer for 



molecular shocks (Cyganowski et al. 2009 2011), 
or an ionized jet ( Araya et al.||2007 ). 

In contrast, in some cases other mechanisms 
such as scattered continuum can also be responsi- 



ble for this emission. Qiu et al. ( 2008 ) attributed 



such emission in a few high-mass star forming re- 
gions to scattered continuum in outflow cavities, 
based on their morphologies and/or excess emis- 



sion at 3.6-/xm. De Buizer and Vacca (2010b have 



made ground-based infrared spectroscopy of two 
Green Fuzzy objects, and conclude that the knotty 
structures in one of the objects are dominated by 
thermal H 2 emission, while the other is associ- 
ated with scattered continuum. Varricatt (2011) 



shows that the bright part of the extended emis- 
sion in the 4.5-/xm and 2.2-/xm continuum position- 
ally match fairly well in IRAS 17527-2439, sug- 
gesting the same origin. He shows that this object 
is also associated with a bipolar jet in H 2 2.122 
/xm emission, but with different distribution from 
the other wavelengths. Smith et al. (20061 ana- 



lyzed their results with spectra obtained by the 
Infrared Space Observatory, and showed that flu- 
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orescent H2 emission also significantly contribute 
to this region in addition to shocked emission in 
the DR 21 outflow. 

Clear understanding of the nature of the emis- 
sion in individual objects would allow us to in- 
vestigate the activity of the associated high-mass 
protostars in detail. While infrared spectroscopy 



is a powerful tool for this purpose ( De Buizer and 
Vacca 2010), it covers significantly limited areas 
of emission regions in many cases. We therefore 
tackle this issue with the observed fluxes, flux ra- 
tios and color-color diagrams based on IRAC data, 
in particular at 3.6, 4.5 and 5.8 /mi. We com- 
pare the colors of Green Fuzzy emission with vari- 
ous cases, including scattered continuum emission 
associated with a low-mass protostar (L 1527 - 
Tobin ct al. 2008), modeled scattered continuum 



emission, shocked emission in low-mass protostel- 



lar jets, thermal H 2 emission calculated by Takami 



et al. (2010), and modeled PAH (Draine and Li 



2007) and fluorescent H2 emission (Draine and 



Bcrtoldi 1996). In Section 2 we describe our tar- 



gets and data. In Section 3 we describe our analy- 
sis and results with flux ratio maps and color-color 
diagrams. In Section 4, we perform comparisons 
between the observed colors and those for the vari- 
ous emission mechanisms described above. In Sec- 
tion 5 we discuss possible contributions from indi- 
vidual mechanisms. In Section 6 we briefly com- 
ment on the emission mechanisms not included in 
this study. In Section 7 we provide conclusions. 

2. Targets and Data 

Our targets are summarized in Table [2] These 
include six nearby (d=2-3 kpc) high-mass proto- 
stars and three bright candidates for high-mass 
protostellar objects reported by I Cygano wski et al.| 
(2008). The selection of the high mass targets 



is based on Bcuther & Shepherd (2005), who 
reviewed previous interferometric observations of 
high-mass protostellar outflows at high angular 
resolutions. These objects are categorized into 
high-mass protostellar objects without any evi- 
dence for ionizing radiation (HMPO), hypercom- 
pact HII regions (HC HIIs) and ultra-compact HII 
regions (UC HIIs). We select two HPMOs (IRAS 
05358+3543, IRAS 16547-4247); three HC H lis 
(G 35.2-0.7 N, G192.16-3.82, IRAS 20126+4104); 
and an UC HII (W 75 N). The additional three ob- 



jects are selected from the catalog of Cyganowski 
et al. (2008). In these objects the emission is fairly 
extended so that we can apply analysis with color- 
color diagrams at different positions (Sections 3 
and 4). These targets are cataloged as "likely 
massive young stellar objects with outflows" by 
Cyganowski et al. (2008), based on the presence 



of extended 4.5 /im emission without confusion 
from nearby point sources and/or problematic sat- 
uration. In addition to the high-mass protostars, 
we also analyze the data for the L 1527 protostellar 
outflow, a well-known low-mass protostar associ- 
ated with scattered continuum in conical bipolar 
cavities (Tobin et al. 20081. Some of the above 



regions show diffuse PAH emission with an exter- 
nal origin. In each of these objects, the spatial 
variation of the PAH emission in the field is much 
smaller than the flux from the object at 3.6, 4.5 
and 5.8 /im, for which we apply the color-color 
analysis in Sections 3 and 4. 

The Spitzer IRAC data in all four bands (3.6, 
4.5, 5.8, and 8.0 /im) were obtained from the 
Spitzer archive. We use the pipeline reduced post- 
BCD (Basic Calibration Data) images provided by 
the archive. The mean FWHMs of the point re- 
sponse functions (PRFs) are 1".66, 1".72, 1".88 
and 1".98, respectively. Figure [I] shows the con- 
ventional three color images: blue, green, red 
for 3.6, 4.5, and 8.0 /im, respectively, as well as 
one-dimensional cut-outs of the intensity distribu- 
tions in all four bands. All the objects are as- 
sociated with Green Fuzzy emission, and the one- 
dimensional profiles show that these have a similar 
shape at 3.6, 4.5, and 5.8 /im. 

Our contrast adjustment for the three colors in 
Figure [T] was made by setting the upper and lower 
limits based on the specified percentage (95, 96, 
97, 98, 99 or 99.5 %, depending on the objects) of 
the flux distribution in the larger areas. As a re- 
sult, stars can appear blue due to the fact that the 
flux is larger at shorter wavelengths; PAH emission 
can appear red due to excessive emission at 8.0 
/im; the remaining extended emission can appear 
green due to relatively low fluxes of stars and PAH 



at this wavelength. While other authors ( Noriega- 



Crespo et aL"l20 04; Marston et al. 2004; Rathborne 
et al.||2005| |Smith et al.|[2006| [Araya et al.||2007| 
Kumar and Grave|2007| [Shepherd et al.|2007| |Ta 



bin et al.|2007l [20081 IQiu et al.|2008||Teixeira et al 



2008 Cyganowski et al. 2008 2009, 2011 Cham- 
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bers et al. |2009j Morales et al. 2009; Zhang and 



Wang 2009| Simpson et al. 2009 Varricatt 2011) 



do not clearly state their color adjustments, we as- 
sume that most, if not all, of them made it in a 
more or less similar manner. 

Despite the "green" color in the three color im- 
age, none of the sources clearly show excess emis- 
sion at 4.5 /im in the one-dimensional cut-outs in 
Figure [TJ Indeed, Chambers et al. 



(20091, who 



coined the term Green Fuzzies, define a quantita- 
tive criteria for their detection as a 4.5 /im to 3.6 
/im ratio of > 1.8 and a 5.8 /im to 4.5 /im ratio 
of < 2.5. Plotting this criteria, one can easily find 
that those with a marginal deficit at 4.5 /jm can 
also be categorized as Green Fuzzies. Through- 
out, the green color in three-color images may not 
always imply the presence of the 4.5-/im emission 
enhanced over the other IRAC bands as stated in 



the literature (e.g., Rathborne et al. 2005 Araya 
eTaTp)07l [De Buizer and Vacca||2010| ~ 

Figure [TJ also indicates the regions in each ob- 
ject where we apply our analysis with flux ratio 
maps and the pixel by pixel color-color diagrams. 
The criterion is described in Section 3. Through- 
out the paper we focus the analysis on the bright- 
est regions (contrast ratio to the brightest pixels 
to be 1/10-1/100) where we expect high signal- 
to-noise. This means that our analysis misses the 
faint outer regions discussed in previous literature 
in the W 75 N dQiu et al.|[2008|) and G 11.92 -0.61 
regions ( |Cyganowski et al.||2008[ |2009| |201l| P 



Imaging observations of H2 at 2.12 /im have 
been made near these bright regions for all the 
nearby high-mass protostar sample. These are:- 



IRAS 05358+3543 (Varricatt et al. 2010), IRAS 



16547-4247 ([Brooks et al.||2003[), G 35.2-0.7 N 
( |Froebrich et al.|2011| |Lee et al ||2012), G 192.16 



3.82 (Indebetouw et al. 2003 Varricatt et al 



2010), IRAS 20126+4104 (Shepherd et al. 2000 



Varricatt et al. 20101), and W 75 N (Davis et al 



1998 Shepherd et al. 2003}. The H 2 2.12 /im emis- 



sion in IRAS 05358+3543, G 192.16-3.82, and W 
75 N are identified as knots, while that in IRAS 
20126+4104 exhibits a collimated jet. These con- 
trast to the fuzzy morphology of the green emis- 



1 Cyganowski et al. ( 2008 1 measured the total flux of G 
11.92-0.61 of 0.334 Jy in the 4.5-fim band. In the area 
shown in Figure [TJ we measure the total flux of 0.257 Jy, 
implying that 77 % of the flux is included. 



sion in Figure [TJ Figure [2] shows the locations 
of the H2 emission in the three-color images. At 
these positions only faint 4.5-/im emission compo- 
nents are seen at some of these positions. The 
H 2 2.12 /im emission in IRAS 16547-4247 and G 
35.2-0.7 N show more complicated structures in 
shocked gas (Brooks et al. 2003; Froebrich et al. 
2011, Lee et al. 2012). These have a significantly 
larger spatial extension (by a factor of 2-3) than 
the bright IRAC emission shown in Figure [TJ and 
most of the flux is distributed outside the region 
indicated. 

Most of the above literature also provide the 
continuum images at ^2 /im, either K(2.2[im), 
K' (2.12 /im), K s (2.15 /im) or 2.14 /im. Their 
distributions in IRAS 05358+3543, G 35.2-0.7 N, 
G 192.16-3.82, IRAS 20126+4104, and W 75 N 
are remarkably different from H2 2.12 /im, but 
to some extent similar to the IRAC emission in 
Figure [TJ In particular, a bipolar structure is ob- 
served in G 35.2-0.7 N and G 192.16-3.82 in both 



IRAC and the continuum at 2-fxm (see Froebrich 
[ et al.|2011[|Indebetouw et al.|2003HVarricatt et al 
|2010[ ). In the west of G 192.16-3.82 both IRAC 
emission and 2-/im continuum show a horn-like 
structure at similar angular scales. In W 75 N, 
both IRAC emission and 2-/im continuum are ex- 
tended to the south-east side of the protostar, in 
contrast to the orientation of the molecular out- 
flow (Shepherd et al. 2003) and the distribution 



of H2 emission in different directions ( Davis et al 
1998[ |Shepherd et al.||2003[ ). The detection of the 



2-/im continuum at these positions indicates that 
the different distribution seen in the bright parts 
of Green Fuzzy emission and H 2 2.12 /im may not 
be attributable to extinction ( Lee et aT[|2012 ). 



Molecular outflows have been observed in most 
of our nearby high-mass protostar sample. G 
35.2-0.7 N and G 192.16-3.82 are associated with 
molecular outflows with similar morphology to the 
direction of elongation in the Green Fuzzy emis- 
sion ( |Gibb et al.||2003[ |Shepherd et al.||1998[ ), but 



the position angles between the Green Fuzzies and 
molecular outflows are different by ~ 30°. In 
IRAS 05358+3543 and W 75 N, the distribution 
of Green Fuzzy emission is remarkably different 
from the molecular outflows observed at millimc- 



ter wavelengths by Beuther et al. 


(2002); Davis 


et al. 


(1998); Shepherd et al. (2003 


) (see also Qiu 


et al. 


2008). For IRAS 16547-4247, Brooks et al. 
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( 2003 ) shows the presence of extended shocked H 2 



emission to the north and south of the protostar. 
This contrasts to the distribution of Green Fuzzy 
emission in Figu re |Tj In IRAS 20126+4104 |Cesa-| 
roni et al. ( 1999 ) shows the presence of a molecular 
outflow in the millimeter HCO + and SiO emis- 
sion with the same orientation as the Green Fuzzy 
emission in Figure [T] but with a slightly smaller 
angular scale. This object is also associated with a 
more extended outflow in the millimeter CO emis- 



sion with a remarkably different orientation ( Shep- 
herd et al.|[2000| . 



3. Flux Ratio Maps and Color-Color Dia- 
grams 

In most of the regions the emission at 8.0 /xm is 
severely contaminated by bright diffuse extended 
emission not directly associated with the proto- 
star. Therefore, we will use the data for 3.6, 4.5, 
and 5.8 /xm, excluding 8.0 /im. To produce accu- 
rate flux ratio maps and color-color diagrams, the 
individual images need to be matched in resolu- 
tion and carefully background subtracted. We first 
convolve each image with the appropriate PRFs to 
try to cancel out the effect of different PRFs at 5.8 
/im from 3.6 and 4.5 /im (note that the PRFs at 
3.6 and 4.5 /im are almost identical — see IRAC 
Data Handbook 3.0). The 3.6 and 4.5 /im images 
are convolved with the PRF at 5.8 /im, and the 5.8 
/xm images with the PRF at 3.6 /im. This yields 
an effective angular resolution of ~3". We then 
subtract the diffuse background emission by mea- 
suring it in the x- and //-directions, and fitting it 
using a linear function. Figure [3] shows how this 
process works for one of our sample. This fitting 
process simultaneously allows us to measure the 
uncertainty of flux, which is dominated by non- 
uniform distribution of diffuse emission. These are 
tabulated in Tabled 

To investigate the origin of the Green Fuzzy 
emission, we use pixel by pixel color-color dia- 
grams with [3.6]-[4.5] and [3.6]-[5.8]. As described 
in Section 2, we focus the analysis on the brightest 
regions where we expect the high signal-to-noise to 
clearly show the tendencies in color-color diagrams 
described later. The criterion for selection is sum- 
marized in Table 2. The selection is made based 
on the 4. 5- /mi flux: >15-<7 for most of the objects. 
For W 75 N and G 298.26+0.74 a larger limit is 



applied to exclude the emission associated with 
other object(s) nearby. For IRAS 20126+4104 we 
apply an additional flux limit (>5-c) for 5.8-/im. 
This allows us to remove artifacts in the color- 
color diagrams that are due to a high uncertainty 
in the background subtraction at this wavelength 
for this object. For L 1527 a flux limit of >25-cr is 
applied to clearly show the linear correlation be- 
tween [3.6]- [4.5] and [3.6]- [5.8] described later. In 
all the objects we rebin the data onto a 3" grid, 
i.e., the same angular scale as the resolution after 
convolution. 

Figure [| shows the 3.6/4.5- and 3.6/5.8- 
/im flux ratios maps in the regions described 
above. This figure shows that the the distribu- 
tions of the 3.6/4.5- and 3.6/5.8-/im flux ratios 
are similar for each object except G 192.16-3.82. 
IRAS 20126+4104, G 35.2-0.7N, G192.16-3.82, 
and L 1527 are associated with stars near the 
boundary. These appear blue in three-color im- 
ages in Figure [l] and exhibit high 3.6/4.5- and 
3.6/5.8-/im flux ratios ratios in the flux ratio 
maps in Figure |4j The color-color diagrams for 
these objects are made excluding these stars. The 
regions after this removal are indicated in Figure 
|4] as well as Figure [IJ In IRAS 20126+4104, the 
3.6/4.5-/im flux ratio map shows another point 
source with high flux ratios in the west of the 
region. This source is absent or marginal in the 
3.6/5.8-/xm map. The inferred color of this source 
and the rest of the region is discussed later. 

In Figure [4] we also plot the positions of the 
protostar (HPMO, HC/UC H II) measured using 
millimeter interferometry. Although the position 
is in the mask of low signal-to-noise or a satu- 
rated region for a few objects, this approximately 
matches the region with the lowest 3.6/4.5- and 
3.6/5.8-/im flux ratios. Similarly, these flux ratios 
are the lowest at the center of L 1527 where the 
protostar is located ( Tobin et al.||2008 ). 

Figure [5] shows the colors measured in the 3x3 
arcsec image bins for Green Fuzzy emission in each 
high-mass star forming region. All the objects but 
IRAS 20126+4104 and G 192.16-3.82 show a fairly 
strong correlation between the [3.6]-[4.5] and [3.6]- 
[5.8] colors along the direction of the extinction 
vector. Their values vary between objects in the 
ranges of [3.6]-[4.5] and [3.6]-[5.8] of 0.4-3 and 0.1- 
4.7, respectively. A similar trend is also observed 
in most of pixels selected in Green Fuzzy emission 
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associated with IRAS 20126+4104. In addition to 
the component in which the [3.6]-[4.5] and [3.6]- 
[5.8] colors show a linear correlation, this object 
is associated with several pixels for which [3.6]- 
[5.8] is constant (2.2-2.4) over [3.6]-[4.5] =0.8-1.5. 
Some pixels in G 192.16-3.82 also show a linear 
correlation like that described above, however, this 
object is also associated with a number of pixels 
with larger [3.6]-[4.5] and [3.6]-[5.8] colors, rang- 
ing between 0.2-1.7 and 1.0-2.4 for [3.6]-[4.5] and 
[3.6]-[5.8], respectively. 

Figure [6] shows the locations of the regions 
which have a [3.6]-[5.8] excess from the linear cor- 
relation in G 192.16-3.82, and where the [3.6]-[5.8] 
color has a constant value in IRAS 20126+4104. 
In G 192.16-3.82 the [3.6]-[5.8] excess is located 
in the outer regions where the 4.5-/mi flux is rel- 
atively faint. In IRAS 20126+4104 the emission 
with a constant [3.6]-[5.8] is located around a point 
source to the west of the protostar. The latter is 
also clearly seen in the 3.6/4.5-/im flux ratio map 
in Figure[3] These areas are responsible for a small 
fraction of the 4.5-/^m flux as compared with that 
of the entire area shown with dashed curves: these 
are 22 % and 6 % for G 192.16-3.82 and IRAS 
20126+4104, respectively. 

Some authors mention the possibility that such 
"green emission" is due to CO ro- vibrational tran- 
sitions at 4.5-5 /Ltm (see e.g., Shepherd et al.|2007" 



Cyganowski et al.| [2008} |2009[ |Chambers et al.| 
2009| , ; see also Takami et al. 2010 and ref- 
erence therein). However, our results do not 
show clear evidence for such emission components, 
which should only appear in the 4.5-/im band. If 
the Green Fuzzy emission in these objects is as- 
sociated with shocked emission as discussed by 
some authors (see Section 1), the absence of CO 
emission contrasts to the shocks associated with 
low-mass protostars, which also show green colors 



in three-color images (e.g., Noriega- Crespo et al. 


2004 


Tobin et al. 2007; Teixeira et al. 2008; Zhang 


and Wang 


20091. In the latter a 4.5-/zm excess, 



presumably due to CO emission, is observed at 
some positions in most of the objects in Takami 
et al. (2010, 2011). In the next two sections we 
discuss the emission mechanism of Green Fuzzy 
emission associated with high- mass protostars and 
their candidates. 



4. Comparisons with Other Objects and 
Models in Color-Color Diagrams 

In the following subsections we compare these 
colors with scattered continuum (Section 4.1), 
molecular shocks associated with low-mass proto- 
stars and modeled thermal H2 emission (Section 
4.2), and PAH and fluorescent H 2 emission (Sec- 
tion 4.3). 

4.1. Scattered continuum 

Figure [7] shows the same color-color diagrams 
as Figure [5] but we include the points for scat- 
tered continuum from the L 1527 protostellar out- 
flow in each diagram. The correlation of [3. 6] -[4. 5] 
and [3.5]-[5.8] colors observed in most of the high- 
mass star forming regions is remarkably similar to 
that of L 1527. The regression line for L 1527 fits 
the observed colors in IRAS 05358+3543, IRAS 
16547-4247, G 298.26+0.74, and most of the col- 
ors measured in IRAS 20126+4104 as well. The 
figure shows a slightly larger [3.6]- [5.8] color (up 
to 0.6) in W 75 N and a slightly smaller [3.6]- 
[5.8] color (by 0.3-0.7) in G 11.92-0.61 and G 
324.72+0.34 compared with the regression line. G 
35.2-0.7 N and W 75 N show a [3.6]-[5.8] excess 
at [3.6]-[4.5]=1.0-1.5. 

If the Green Fuzzy emission in our sample of 
high-mass star forming regions is due to scat- 
tered continuum, the observed colors should be 
a function of: (1) the intrinsic color of the central 
star(+disk); (2) the color change due to scatter- 
ing; and (3) extinction. To investigate the first two 
issues, we have made simplified calculations with 
existing star-disk models and dust models. The 
modeled colors for star-disk systems were obtained 
from Robitaille et al. (20061, based on radiative 



transfer calculations for protostellar environments 
with the protostar, disk, envelope and outflow cav- 



ities by Whitney et al. (2003b). Robitaille et al 



(2006) present a grid of infrared spectral energy 
distributions (SEDs) for 200,000 cases, with a va- 
riety of stellar and disk parameters, at ten viewing 
angles for each model. To derive the infrared flux 
from the disk, these authors include "active vis- 
cous heating based on the standard accretion disk 
model by Shakura and Sunyaev ( 1973 ) , and "pas- 



sive heating by stellar radiation. They also adopt 



the disk geometry of Shakura and Sunyaev ( 1973 1 



To derive the intrinsic color of the star-disk sys- 
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tems (i.e., "1" described above), we select all the 
modeled results with Ay < 5 (corresponding to 
^3.6/jm-4.5/jm < 0.07 and -E3.6yum-5.8/zm < 0.12 
for the adopted grain model). This criterion is 



applied due to the fact that the Robitaille et al 



(2006) models apparently include dust in the en- 



velope and outflow cavity so that we cannot derive 
the values without extinction. The above criterion 
allows us to include a large number of samples 
with a variety of disk masses and accretion rate 
with small errors in color due to extinction. We 
also select all the results with a stellar mass larger 
than 8 M G . As a result, a total of 17,580 SEDs 
are selected from the grid, with a range of stellar 
mass M*=8-50 M , stellar effective temperature 
T ~ 4 x 10 3 - 5 x 10 4 K, disk mass 0-3 M , and 
accretion rate M ~ 10~ 14 - 10~ 3 M yr _1 if the 
system hosts a disk. In addition to the modeled 



colors of Robitaille et al. (2006), we calculate the 



color for blackbodies with a single temperature of 
T=500, 1000 and 4000 K. 

The color change via single scattering is derived 
from Robitaille et al. ( 2006 1 , with the grain models 



based on the optical constants and size distribu- 



tion by Laor and Draine (1993) and Kim et al 



(1994), respectively. We consider the following 
two cases for scattering: (1) scattering in opti- 
cally thin regions, i.e., where the color change via 
scattering is determined by the scattering cross 
section at each wavelength; and (2) scattering for 
optically thick and geometrically thin cavity walls, 
which may be realistic for outflow cavities associ- 



ated with low-mass protostars (e.g., Whitney et al. 
2003al iTobin et al.||2008|). For the latter, the effi- 



ciency of the scattering should be approximately 
determined by the scattering albedo, assuming 
that the effects of internal extinction and multi- 
ple scattering are negligible. For the color change 
due to extinction, the extinction vector based on 



Chapman et al. ( 2009 1 has a large uncertainty as 
shown in Figure [5j Here we assume that the linear 
correlation observed in L 1527 is due to different 
extinction at different positions, and extract the 
values from the regression line. 

Figure [8] shows the modeled colors of the star- 
disk systems with scattering together with the ob- 
served colors. This figure shows that colors of 
these objects approximately lie within the extinc- 
tion track of either the optically thin or thick 
results. In this case, the linear correlation of 



the [3.6]-[4.5] and [3.6]-[5.8] colors in each ob- 
ject is attributed to different extinctions at dif- 
ferent positions. The different objects show dif- 
ferent colors not only in the direction of the ex- 
tinction vector, but also across it. The latter is 
consistent with the idea that the intrinsic color 
of the star(+disk) systems differs between ob- 
jects. In IRAS 05358+3543, G 11.92-0.61, and 
G 324.72+0.34, the intrinsic flux of the central 
source is dominated by a star with a relatively 
high color temperature (T >1000 K). In W 75 
N and G 298.26+0.76, the central source is red- 
der, i.e., its flux is dominated by a protostar with 
a relatively low color temperature (T <1000 K). 
Figure [8] also shows that the two different scatter- 
ing models yield similar results for the color of the 
embedded central source. 

In G 35.2-0.7 and W 75 N the distribution of 
the correlation is broader, the [3.6]-[5.8] color par- 
tially exceeding those expected for the models at 
[3.6]-[4.5] ~1. Figure [9] shows their distribution 
in the 4.5-^tm map. These areas are responsible 
for 34 % and 9 % of the 4.5-/im flux integrated 
over the entire area shown in dashed curves in G 
35.2-0.7 and W 75 N, respectively. 

4.2. Shocks in low-mass protostellar jets, 
thermal H2 emission 

Figure [10| shows the same color-color diagrams 
as Figures [5j and [7] but with the colors observed in 
six low-mass protostellar jets analyzed by Takami 
et al. (2010). They showed that while some of 



them are explained by thermal H 2 emission (a 
combination of rotational and ro- vibrational tran- 
sitions), others requires additional emission at 4.5 
/im, presumably due to the vibrationally excited 
CO emission. In Figure [lOj the colors are dis- 
tributed along an arc with [3.6]-[4.5] and [3.6]- 
[5.8] of 1.1-2.3 and 1.4-3.5, respectively, located 
at both sides of the linear color correlation ob- 
served in most of the objects. While the observed 
correlations between [3.6]-[4.5] and [3.6]-[5.8] are 
totally different between Green Fuzzy emission 
and shocks associated with low-mass protostars, 
their colors match in some color ranges in most 
of the objects. The range of such colors are [3.6]- 
[4.5] = 1.2-1.6 and [3.6]-[5.8] = 1.5-2.5 for IRAS 
05358+3543, G 192.16-3.82, IRAS 20126+4104, 
W 75 N and G298.26+0.74. This color range is 
slightly larger in G 35.2-0.7 N ([3.6]-[4.5] = 1.1-1.6 
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and [3.6]-[5.8] = 1.6-2.8), while IRAS 20126+4104 
has an additional overlap at [3.6]- [4.5] = 1.3 and 
[3.6]-[5.8] = 2.4. 

Figure [TT] shows the colors for thermal H2, i.e., 
the primary source of H2 emission in shocks, based 
on calculations by Takami et al. (2010). As in 



Takami et al. (20101 the models are made for 
iso-thermal cases and shock slabs with a power- 
law cooling function (A oc T~ a ). The latter has 
been used to explain shocked H2 emission ob- 



served both from the ground (e.g., Brand et al. 


1988 |Gredel 


1994 Everett et al.||1995 Richter 


et al. 1995 


Takami et al. 


2006) and in space 


(Neufeld and Yuan 


2008 


Neufeld et al. 2009 


Takami et al. 2010 


For calculations of non- 



local-thermal-equilibrium (non-LTE) cases ther- 
mal collisions are included for two cases, H+He 
and H2+He, corresponding to gas with relatively 
high (> 0.002 - .02) and low (< 002 - 0.02) 
dissociation rates ( Takami et al.||2010 ). We adopt 



A-coeffcients provided by |Wolniewicz et al. ( 1998J), 



and collisional rate coefficients for H2 and He by 



Le Bourlot et al. (1999). For collisional rate co- 



efficients of H, we adopt Wrathmall et al. (2007) 



and Le Bourlot et al. (1999), and show the results 



separately. According to Wrathmall et al. ( 2007 1 , 



they provide rate coefficients with a better accu- 



racy than Le Bourlot et al. ( 1999 1 due to improved 



representation of the vibration eigen functions. In 
contrast, the coefficients provided by |Le Bourlot] 



et al. (1999) can explain the observed IRAC col- 



ors better in Takami et al. (2010). 



The energy levels we include for calculations are 
245 {E u /k = - 430 00 K) for LTE based on al l 
the levels included in Draine and Bertoldi (1996); 



the lowest 49 energy levels (46 lines in the IRAC 
four bands, E/k up to 2 x 10 4 K) for non-LTE 
H2 with H+He collisions; and the lowest 36 en- 
ergy levels (32 lines in the IRAC four bands, E/k 
up to 1.69 x 10 4 K) for non-LTE H 2 with H 2 +Hc 
collisions. The number of transitions included for 
non-LTE calculations is limited by the availability 
of collisional rate coefficients. The level popula- 
tions of ortho- and para-H2 are calculated sepa- 
rately, and those fluxes are combined assuming an 
ortho/para ratio of 3. The upper limit of tempera- 
ture is set at 4000 K, approximately corresponding 
to the dissociation temperature of H2 molecules 
(e.g., Lepp and McCray||1983 ). We set the lowest 
temperature of 30 K for numerical integration of 



the temperature slabs for power-law cooling. Since 
the emission in IRAC bands should originate from 



gas at much higher temperature (>1000 K |Takami 
et al. 2010), this limit does not affect the colors 



we show in Figure 11 For all the above cases the 
measured spectral response functions are used to 
obtain the IRAC fluxes. 

Due to the limited number of transitions in- 
cluded, the accuracy of non-LTE calculations are 
low for the IRAC 3.6-/im band at the highest tem- 
peratures. The gaps between such calculations 
and those with 245 levels are also shown in Fig- 
ure [TT] for the LTE regimes. The non-LTE gas 
of our calculations does not show smaller [3.6]- 
[4.5] and [3.6]-[5.8] colors than LTE, since the 3.6- 
/j,m emission requires a higher density for thermal- 
ization than the emission in the other IRAC bands 



dTakami et al.||2010[ ). 

In Figure [12] we plot the observed colors in 
Green Fuzzy emission with possible colors due to 



thermal H2 as shown in Figure 11 While thermal 
H2 could explain the colors observed at most of 
the positions, it cannot explain the colors with 
small [3.6]-[4.5] and [3.6]-[5.8] colors (<1.0 and 
<1.5, respectively) observed in IRAS 05358+3543, 
IRAS 16547-4247, G 35.2-0.7 N, G 192.16-3.82, 
W 75 N and G 298.26+0.74. These correspond 
to the regions with lowest 3.6/4.5- and 3.6/5.8- 
/im flux ratios in Figure [4] Furthermore, thermal 
H2 emission with H+He collisions cannot account 
for the large [3.6]-[4.5] and [3.6]-[5.8] colors (>2.1 
and >2.6, respectively) in IRAS 16547-4247, W 75 
N, G 11.92-0.61 and G 324.72+0.34. In Figure[TJ 
we compare the linear correlation of the observed 
[3.6]-[4.5] and [3.6]-[5.8] colors with modeled ther- 
mal H2 emission for collisions with H2+He (Figure 
11 1 in more detail. The figure shows that the linear 
correlation observed in Green Fuzzy emission may 
be attributed to different densities of gas observed 
in thermal H 2 emission. The problems with this 
explanation are discussed in Section 5 in detail. 

4.3. PAH and fluorescent H2 emission 

Figure [14] shows color-color diagrams for the 
models of PAHs and fluorescent H2 emission, 
usually associated with photodissociation regions 
(PDRs; 



e.g. 



Hollcnbach and Tielens 1997 for a 



review). Their fluxes are obtained from Draine 
and Li| ( |2007[ ) and |Draine and Bertoldi| (fl996|) 
respectively. The models by Draine and Li (2007) 
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also include thermal dust continuum, which is 
prominently observed at longer wavelengths, but 
its contribution is negligible in the IRAC bands as 
compared with PAHs. For comparison we selected 
their results for a UV field 10 2 -10 5 as large as the 
interstellar radiation field for the solar neighbor- 
hood. Such a UV field is comparable to those 
observed in well-known dense PDRs associated 
with OB stars such as Orion Bar, NGC 2023 and 



7023 (e.g., 


Tielens 2008 


Usuda et al. 1996 


Takami 


et al. |2000). For the abundance of PAH, all the 



results calculated for the dust/gas ratio in the 
Milky Way (fraction of C of 0.47-4.58 % to the 
entire abundance in the interstellar medium) are 
included in the figure. 



In Figure 14 the colors for PAHs are noticeably 
offset from the linear correlation of [3.6]-[4.5] and 
[3.6]-[5.8] colors, with a combination of small [3.6]- 
[4.5] colors and large [3.6]-[5.8]. This indicates 
that PAHs are not primarily responsible for the 
color distribution observed in Green Fuzzy emis- 
sion. In contrast, the constant [3.4]-[5.8] colors 
observed in several positions in IRAS 20126+4104 
are similar to those of PAHs. The figure shows 
that the modeled colors for fluorescent H2 show 
similar correlation between the [3.6]-[4.5] and [3.6]- 
[5.8] colors to those observed in most of Green 
Fuzzies. The discrepancy between the modeled 
and observed colors could be attributed to an ex- 
tinction Ax up to ^20. This combination, how- 
ever, cannot explain some pixels in G 35.2-0.7 N 
and G 324.72+0.34, with the lowest [3.6]-[4.5] and 
[3.6]-[5.8] colors. We discuss the possible contri- 
bution of fluorescent H2 emission to the observed 
colors in detail in Section 5.3. 

5. Discussion 

5.1. Scattered continuum as a possible pri- 
mary source of Green Fuzzy emission 

The colors of the star+disk flux via scatter- 
ing and extinction are approximately consistent 
with the observed colors. In particular, a sim- 
ilar correlation between [3.6]-[4.5] and [3.6]-[5.8] 
colors are observed in scattered continuum in the 
cavity of the L 1527 protostellar outflow, agree- 
ing with this explanation. In this case, the color 
distribution observed in each object can be at- 
tributed to differing foreground extinctions (in- 
cluding dense cores and envelopes). This is cor- 



roborated by the facts that (1) the position of 
the protostar matches the regions with the low- 
est I 3 . 6 [im/h.5 ^m and I 3 . 6 ^m/h.s fim ratios 
for most of the objects (Figure |4J); and (2) the 
column density of the gas and dust is higher 
close to the protostar than in surrounding re- 
gions (see references for Table |2j. This also 
suggests that the internal extinction in the ex- 
tended emission regions is relatively small, other- 
wise the J3.6 fim/h.5 fim and I 3 , 6 fim/h.8 fim ra- 
tios would decrease with increasing distance from 
the protostar due to reddening. In this context, 
the Green Fuzzy emission we analyzed may repre- 
sent the morphology of cavities due to outflowing 
gas, and the J 3 . 6 fj,m/h.5 fim and J 3 . 6 fim/h.8 fim 
ratio maps in Figure [4] may represent the distribu- 
tion of foreground extinction. 

The explanation that Green Fuzzy emission is 
scattered continuum associated with outflow cavi- 
ties is consistent with the fact their distribution is 
more similar to the 2-/im continuum than to the 
H2 2.12 fim emission (Section 2). The morpho- 
logical discrepancies between the 4.5-^tm emission 
and continuum at 2 fim could be attributed to dif- 
ferent extinctions at these wavelengths. This ex- 
planation is also consistent with the fact that such 
emission in some objects are associated with signa- 
tures of molecular outflows and/or outflow shocks 
(Section 2; see also Shepherd et al. 2007 Araya 



et al.|2007trQiu et al. 



2008; Cyganowski et al.|2009| 



2011). It is intriguing that the distribution of the 
IRAC emission and/or 2- fim continuum do not 
match H2 emission and/or molecular outflows in 
some objects, including IRAS 05358+3543, IRAS 
16547-4247 and W 75 N (Section 2). This would 
indicate the complicated nature of high-mass star 
forming regions, which are often associated with 
multiple protostars and outflows like these objects 



(e-g-, 


Shepherd et al. 


2000 


2003 


Beuther et al. 


2002 


Franco-Hernandez et al.||200£ 


)■ 



As shown in Sections 3 and 4.2, some positions 
in G 192.16-3.82, G 35.2-0.7 N, and W 75 N show 
a [3.5]- [5.8] excess from the linear correlation of 
the [3.6]-[4.5] and [3.6]-[5.8] colors in the remain- 
ing region. This may be attributed to a different 
origin than scattered continuum. The deviation of 
the color observed in G 192.16-3.82 approximately 



matches the colors of PAHs in Figure 14 suggest- 
ing that this is due to diffuse contaminating PAH 
emission in the outer region where Green Fuzzy 
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emission is relatively faint (Figure [6]). This expla- 
nation is corroborated by the fact that the objects 
near G 192.16-3.82 are associated with extended 
emission at 8.0-/xm (Figures 1 and 2), a signpost 
of PAH emission (Section 1). The same explana- 
tion may also apply to the [3.6]- [5.8] excess in W 
75 N. This is corroborated by the fact that the re- 
gion associated with the [3.6]-[5.8] excess is close 
to extended 8.0-/xm emission (Figures 1 and 2) as- 
sociated with an H II region ( Haschick et al.|1981[ ). 

An alternative possible contributor for W 75 N 
and G 35.2-0.7 N is shocks the physical conditions 
of which are similar to some of those associated 
with low-mass protostars. This is consistent with 
the fact that the [3.6]-[5.8] excess in these objects 
is observed over color ranges similar to such shocks 
in Figure 
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Note that the fraction of the 4.5- 
fim flux in this component to that of the entire re- 
gion (22, 34, 9 % for G 192.16-3.82, G 35.2-0.7 N, 
and W 75 N; Sections 3 and 4.2) is an upper limit 
for the additional emission component, assuming 
that the flux in the marked regions is solely due to 
shocks or PAHs. Throughout, the contaminating 
emission is a minor component in each object in 
terms of both spatial coverage and flux. 

5.2. Thermal H2 emission as a possible 
primary source of Green Fuzzy emis- 
sion 

As described in Section 1, H2 or CO emis- 
sion in shocks is often regarded as the primary 



source of Green Fuzzy emission ( Rathborne et al. 
2005} [Aray^"eTain[20C]7l |Shepherd et al.|2007||Qiu' 
et al.|2008HMorales et al.||2009"]|Cyganowski et al 
2008H2009n2011| IChambers et al.||2009| ). Such 
emission is usually expected via thermal excita- 
tion. Indeed, thermal H2 and CO emission explain 
well the IRAC colors observed in shocks associ- 



ated with low-mass protostars (Ybarra and Lada 
20091 i Takami et al.|2010] ). However, the [3.6]-[4.5] 



and [3.6]-[5.8] colors show different correlations be- 
tween Green Fuzzy emission and these shocks asso- 
ciated with low-mass protostars (Section 4.1, Fig- 
ure 10). This implies that different shock condi- 



tions are required if Green Fuzzy emission is due 
to shocks. 

As described in Section 4, the Green Fuzzy 
emission in our sample does not show any clear 
signature for the presence of the ro-vibrational 
CO emission, which should only appear in the 4.5- 



/im band. In this subsection we therefore focus our 
discussion on thermal H2 emission, which is pri- 
marily responsible for shocked H 2 emission in most 
cases (e.g.,|Brand et al.|1988||Gredel|1994||Everett 
[et~al~1[T995l |Richter et al.|1995||Usuda et al.||1996 



Eisl offel et al.|2000| |Takami et al.|2006||Beck et aT 



20081 |Neufeld et al. 2006[ |Neufeld and Yuan|2008 



Neufeld et al. 2009). Indeed, the colors observed 
at most of the positions in Green Fuzzy emission 
are consistent with thermal H2 emission (Section 



4.1, Figure 12), which is primarily responsible for 
shocked H2 emission. Furthermore, the linear cor- 
relation of the [3.6]-[4.5] and [3.6]-[5.8] colors may 
be attributed to different densities if thermal H 2 
emission with H2 and He collisions is responsible 
for the Green Fuzzy emission (Section 3.1, Figure 
131. 

These explanations, however, face several prob- 
lems described below. First, in six out of nine 
objects (IRAS 05358+3543, IRAS 16547-4247, G 
35.2-0.7 N, G192.16-3.82,W 75 N, G298.26-0.74), 
some positions show smaller [3. 6] -[4. 5] and [3.6]- 
[5.8] colors (<1.0 and <1.5, respectively) than 
those predicted by thermal H 2 emission (Section 
4.1, Figure |l2"|) . In these objects the linear cor- 
relation between [3. 6] -[4. 5] and [3.6]- [5.8] is ob- 
served across both color ranges, suggestive of an- 
other origin. Secondly, Figures [12] and [13] would 
indicate the temperature of thermal H 2 is ^2000 
K, in which H 2 should be efficiently excited for 



the 1-0 S(l) 2.12 /im emission (e.g 


, Brand et al. 


19881 |Gredel| 


|1994t |Everett et al.| 


|1995| |Richter 


et al.|1995[|Usuda et al.|1996t|Eisloffel et al. 2000; 


Takami et al.|2006 Beck et al.|2008 


Neufeld et al.] 


2006 Neufeld and Yuan 2008 


Neufeld et al. 2009). 



It would therefore be puzzling if the distribution 
of Green Fuzzy emission were remarkably different 
from H2 2.12 /im emission (Section 2). See also 



Takami et al. (2010) for the dependence on den- 
sity and temperature of H 2 emission at the 2.12 
fj,m and IRAC bands. 

Thirdly, if this color correlation is due to the 
different densities in Figure [13] it implies that 
the positions with larger [3.6]-[4.5] and [3.6]- 
[5.8] colors (thereby lower /3.6/^m/^4.5/zm an d 
I3.6fj.rn/ hs/im flux ratios) are associated with 
lower densities. Since the /3.6/im/^4.5/im and 
h.Gfxm/ hs/im flux ratios increases with distance 
in G 35.2-0.7 N, G 192.16-3.82, W 75 N, the 
density should increase with distance in these ob- 
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jects. This is in the opposite sense as one would 
expect in outflow shocks. If these shocks are due 
to interaction between the ejecta and surrounding 



gas, as in many case (e.g., Bally et aL||2007 Arce 



et al. 2007), one would expect lower densities in 



the downstream as for the distribution of molecu- 



lar gas around protostars (e.g., Stahler and Paha 
2005). If these shocks are in the ejecta, one would 



also expect such a density distribution as long as 
the flow has an opening angle similar to that ob- 



served in many molecular outflows (e.g., Beuther 
and Shepherd|[2005{ |Arce et al.|2007 >. 



Fourthly, it is likely that the ro-vibrational CO 
emission significantly contributes to the entire 4.5- 
/im flux if the hydrogen number density exceeds 
~ 10 7 cm" 3 (iNeufeld and Yuan|2008[|Takami et al 



2010). This would cause larger [3.6]-[4.5] colors 



at high densities, as observed in colors in shocks 
associated with low-mass protostars at [3.6]-[5.8] 

for 



sec; 



Takami et al. 2010 



-1.5-2.0 (Figure [10 
details). However, such a tendency is not seen for 
the colors of Green Fuzzy emission in Figure |13[ 

Throughout, we conclude that it is not likely 
that Green Fuzzy emission is primarily due to 
thermal H2 emission. One might think that a 
combination of thermal H2 emission and scattered 
continuum could explain the linear correlation we 
have discussed. This explanation could overcome 
the first problem described above, interpreting the 
lowest [3.6]-[4.5] and [3.6]-[5.8] colors as scattered 
continuum. However, the second problem still 
remains. The third and fourth problems could 
only be overcome if the shock conditions were uni- 
form, yielding a single combination of [3.6]-[4.5] 
and [3.6]-[5.8] colors over the region. This is not 
likely, considering the fact that shocks associated 
with high-mass protostellar outflows show compli- 
cated shock structures (see, e.g. 



Kaifu ct al. 2000 



Davis et al. 2007 Cunningham et al. 2009), sug- 



gesting the presence of different physical condi- 
tions in each outflow. It is noteworthy that, ac- 
cording to this explanation, we would expect scat- 
tered continuum further away from the protostar 
than H2 emission (i.e., the regions with the largest 
^3.6/xm/-^4.5//m an d ^3.6/im/-fs.8/Ltm) m G 35.2-0.7 
N. This is the opposite trend as observed in the 
H 2 2.12-^m emission and scattered continuum at 
2 /zm (iFroebrich et al.|2011l|Lee et al.|2012|). There 



is no clear theory for how the above explanation 
overcomes this discrepancy. 



Alternatively, a combination of thermal and flu- 
orescent H2 emission in shocks may explain the 
linear correlation in the color-color diagrams. We 
discuss this in the next subsection. 

5.3. Fluorescent H2 

The linear correlation of the [3. 6] -[4. 5] and [3.6]- 
[5.8] colors may also be explained by a combi- 
nation of fluorescent H2 emission and extinction 
(Section 4.3). To investigate this possibility in de- 
tail, we calculate the IRAC fluxes expected for flu- 
orescent H 2 in PDRs based on |Draine and Bertoldi| 
(fT996l. As for Section 4.3 and Figure 
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1996 



Table 



we 

selected their results for a UV field 10 2 -10 5 as 
large as the interstellar radiation field for the solar 
neighborhood. Again, such a UV field is compara- 
ble to those observed in well-known dense PDRs 
associated with OB stars such as Orion Bar, NGC 
2023 and 7023 (e.g., |Tielens||2008| |Usuda et al 
Takami et al.||2000| ). These are tabulated in 
3} The table shows that the predicted fluxes 
are up to —3 MJy str -1 , significantly lower than 
those observed in the bright areas of the Green 
Fuzzy emission shown in Figure [l] (> 10 MJy 
str- 1 ). 

Furthermore, PDRs are also associated with the 
PAH emission, which appears to be much brighter 
in the 3-8 /itm range (Tielens 2008). To inves- 



tigate this issue more quantitatively, we calculate 
the modeled fluxes for the same range of UV fields 



based on Draine and Li (2007). Draine and Li 



(2007) provide fluxes per column density, with- 
out the models for the PDRs. We therefore as- 
sume a column density for the PDRs of Ay=l, 



a typical column dens ity for PDRs (e.g., Hollen- 
bach and Tielens 1997), corresponding to 1.9x 10' 21 
cm~ 2 ( Mat his 2000). These are also tabulated 
111 Table |3| for minimum and maximum abun- 



dances of PAHs modeled by Draine and Li ( 2007 ) 
(qpAH=0A7 and 4.58 %, respectively). The larger 
qpAH value is comparable to those of well-studied 

Table [| 



dense PDRs (~3.5 %, Tielens 2008). 



shows that, even for the minimum abundance of 
PAH (i.e., lower than those of well studied dense 
PDRs by a factor of ~7), the modeled PAH fluxes 
are larger than fluorescent H2 by a factor of >6 
and >100 at 3.6 and 5.8 /xm, respectively. Thus, 
despite the uncertainty in the calculations (i.e., 
the column density of the PDRs), it is not likely 
that the fluorescent H 2 emission dominates over 
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the PAH emissions in normal PDRs. 

Fluorescent H 2 emission has been observed in 



the high-mass protostellar outflow DR 21 (Fernan- 



des et al. 1997 Smith et al. 2006 ) and a well-known 



Herbig-Haro object ( Fernandes and Brand| [l995 ) . 
A combination of thermal and fluorescent H2 may 
consistently explain the linear correlation of [3.6]- 
[4.5] and [3.6]- [5.8] colors observed in Green Fuzzy 
emission in most of the objects (Figures 12] 13] 
14), if most of the PAHs+dust grains are destroyed 



by the passage of shocks. Fernandes et al. (1997) 



and |Smith et al.| ([2006]) suggested a UV field of H 2 
excitation of X ~ 1CP~- 10 4 in the DR 21 outflow. 
Again, such a UV field would yield H 2 flux too 
low compared to the observations. This problem 
may be overcome if (1) most of the dust grains 
are destroyed by shocks, allowing UV photons to 
be absorbed more efficiently by molecular hydro- 
gen; and/or (2) turbulence in the shocks make the 
ultraviolet H2 lines broader, allowing more pho- 
tons to be absorbed. Models for the coexistence 
of fluorescent and thermal H2 in shocks are re- 
quired to further investigate this possibility. Alter- 
natively, follow-up spectroscopic observations of 
Green Fuzzy emission would allow us to further 
constrain the associated shock conditions. 

5.4. PAHs 

Our color-color diagrams show that PAHs are 
not the primary emission mechanism for Green 
Fuzzy emission (Section 4.3, Figure 14). In this 



section we discuss the cases for IRAS 20126+4104, 
in which a relatively constant [3.6]-[5.8] color is ob- 
served at several positions. 

These positions in IRAS 20126+4104 show a 
relatively constant [3.6]- [5.8] color similar to PAHs 
(Section 4.3, Figure 14). While we do not reject 



the possibility that these emission components are 
primarily due to thermal H 2 (Section 4.1, Fig- 
ures 12 13 1, this coincidence suggests that at least 
their 3.6- and 5.8 /im fluxes are primarily due to 
PAHs. However, the observed [3. 6] -[4. 5] colors in 
these regions are significantly larger than those of 
PAHs modeled by Draine and Li ( |2007 ), ranging 
from 0.8 to 1.5. This may be primarily due to a 
contribution from another emission mechanism at 
4.5-/xm, in which the PAHs emission is the faintest 
among the four IRAC bands (e.g., Reach et al 



[20061 [Draine and Li|[2007| |Tielens||2008[ see also 
Table ]3j). 



In IRAS 20126+4104, this emission component 
is associated with a point source to the west of the 
high-mass protostar (Section 3, Figure [6]. The 
three-color image in Figure [2] shows the presence 
of excess emission at 8.0-/xm, supporting the idea 
that this region is associated with PAH emission 
(Section 1). This region is surrounded by emission 
with different colors, which can be attributed to ei- 
ther scattered continuum (Section 5.1) or H 2 emis- 
sion (Sections 5.2, 5.3). Thus, the different [3.6]- 
[4.5] colors described above can be naturally at- 
tributed to different contributions from two emis- 
sion components (i.e., PAH and scattered contin- 
uum, or PAH and H 2 emission). 

Four of our targets are associated with hyper- 
compact and ultracompact H II regions (Table |]|. 
Thus, these should also be associated with far- 
UV radiation which can excite PAHs. However, 
none of our targets show clear evidence for the 
presence of PAH emission directly associated with 
the protostar. These suggest that the high-mass 
protostars and hypercompact/ultracompact H II 
regions are heavily embedded even at the wave- 
lengths of our analysis, i.e., 3.6-5.8 /xm. 

6. Caveats 

It is beyond the scope of this paper to discuss 
the emission mechanisms below, but we briefly 
state their possible contributions. 

Thermal dust continuum might also contribute 
to the 3.6-5.8 /xm flux. Such emission can be ex- 



tended in the outflow cavity (e.g., De Buizer 2006 1 



According to Draine and Li (2007) this emission 



component is negligible compared with PAH emis- 
sion in the IRAC bands. However, the contribu- 
tion from the thermal dust continuum may not be 
negligible in circumstances where PAHs are de- 
stroyed but larger grains survive. 

Fast shocks should be associated with atomic 
and ionic lines, such as Pf 7 3.74 /xm, Br a 4.05 
/xm 



et al. 2006 



Pf P 4.65 /xm, and Fc II 5.34 /xm ( |Reach 
see Table [T] for details), and these 



may also contribute to the 3.6-5.8 /jm bands. In 



the limited cases discussed by Reach et al. ( 2006 ) 



the above lines in shocks yield J3.6/xm/^4.5/xm and 
^3.6/xm/^5.8/xm ratios of < 0.1, significantly lower 
than those observed in Green Fuzzies (Figure [4| . 
Thus, this possibility should be addressed with 
more detailed modeling and/or spectroscopic ob- 
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ser vat ions. 

7. Conclusions 

We analyzed the flux ratios and pixel-pixel col- 
ors of "Green Fuzzy" emission toward six nearby 
(d=2-3 kpc) well-studied high-mass protostars 
and three candidate high-mass protorstars ob- 
served using the archival data for Spitzcr IRAC 
at 3.6, 4.5, and 5.8 /an. In color-color diagrams 
most of the sources show a positive correlation be- 
tween the [3.6]-[4.5] and [3.5]-[5.8] colors, in all or a 
part of the region, along the extinction vector. We 
compare them with the colors of other objects and 
models, i.e., (1) modeled scattered continuum, and 
observed scattered continuum associated with the 
L 1527 low-mass protostellar outflow; (2) shocks 
associated with low-mass protostars and modeled 
thermal H 2 emission; (3) modeled emission for flu- 
orescent H 2 ; and (4) modeled emission for PAHs. 
The conclusions are as follows:- 

1. Of the above emission mechanisms, scat- 
tered continuum in the outflow cavities pro- 
vides the simplest explanation for the ob- 
served linear correlation. Indeed, such a cor- 
relation is also observed in the scattered con- 
tinuum in L 1527. In this case, the differ- 
ent colors within each object are attributed 
to different degrees of foreground extinction 
at different positions in the extended emis- 
sion. Different objects show different colors 
not only along the extinction vector, but also 
across it. This can be attributed to different 
intrinsic colors of the star(+disk) system be- 
tween the objects. This interpretation, that 
scattered continuum is responsible for Green 
Fuzzy emission, is consistent with the fact 
that the distribution of emission is remark- 
ably different from H 2 2.12 /im emission in 
some objects, but more similar to the con- 
tinuum at 2 /jm. 

2. The observed color correlation is remark- 
ably different from that observed in molecu- 
lar shocks in low-mass protostars, in which 
the emission is due to thermal H 2 , plus 
CO at high densities. Even so, a signifi- 
cant fraction of the Green Fuzzy emission 
show colors within the color range of mod- 
eled thermal H 2 emission. However, this 



emission mechanism is not likely for the fol- 
lowing reasons: (1) the observed linear cor- 
relation in some objects exceeds the color 
range predicted for thermal H 2 ; (2) while ro- 
vibrational CO emission should contribute 
to the 4.5-/um emission at high densities 
( n H <^ 10 7 cm ~ 3 ), the observed color-color 
diagrams do not clearly show such a signa- 
ture; (3) the temperature inferred by the 
color diagrams (^2000 K) would allow for 
H 2 2.12^m emission to be observed, how- 
ever, the distribution of this line emission 
does not match Green Fuzzy emission in 
some objects; and (4) higher densities in the 
outer region compared to the inner region 
are required to explain the color correlation 
in some objects, however, this is not likely. 
The first problem may be overcome if we at- 
tribute the lowest [3. 6] -[4. 5] and [3. 6] -[5. 8] 
colors to another origin, e.g., scattered con- 
tinuum. However, a combination of thermal 
H 2 emission and scattered continuum does 
not seem to overcome the other problems de- 
scribed above. 

3. Fluorescent H 2 in normal PDRs cannot ac- 
count for the observed flux. Furthermore, 
emission from PAHs should dominate over 
fluorescent H 2 in the IRAC bands in such 
circumstances. Fluorescent H 2 in shocks 
may overcome these problems if shocks de- 
stroy PAHs and dust grains, and turbulence 
in shocks allows far-UV photons to be ab- 
sorbed by H 2 more efficiently. A combina- 
tion of fluorescent and thermal H 2 may also 
explain the observed linear correlation be- 
tween [3.6]-[4.5] and [3.6]-[5.8] colors. Shock 
models with these excitation mechanisms are 
required to investigate the feasibility of this 
explanation. 

4. As expected, our color-color diagrams show 
that PAHs are not the primary mechanism 
for Green Fuzzy emission. PAHs may sig- 
nificantly contribute in a part of the emis- 
sion regions in another few objects including 
IRAS 20126+4104, G 192.16-3.82 and W 75 
N. None of our sample show clear evidence 
for PAH emission directly associated with 
the high-mass protostars, several of which 
should be associated with ionizing radiation. 
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This suggests that those protostars are heav- 
ily embedded even at mid-infrarcd. 

In summary, among the emission mechanisms 
discussed above, scattered continuum in outflow 
cavities provides the simplest explanation for the 
observed linear correlation. Alternative possible 
emission mechanisms to explain the above linear 
correlation may be a combination of thermal and 
fluorescent H2 emission in shocks, and a combina- 
tion of scattered continuum and thermal H 2 emis- 
sion, but detailed models or spectroscopic follow- 
up are required to further investigate this possi- 
bility. The first interpretation (i.e., attributing 
the Green Fuzzy emission to the scattered con- 
tinuum in outflow cavities) does not change inter- 
pretations of some previous literature which show 
that Green Fuzzy emission is associated with sig- 
natures of jets and outflows (e.g., molecular out- 
flows, shock tracer, or collimated ionized jet). 

This work also highlights the difficulty in ob- 
tain definitive conclusions for the excitation mech- 
anism of Green Fuzzies. Spectroscopic follow-up 
would allow us to complement this study in order 
to understand the nature of high-mass protostars 
themselves, or conditions of outflow shocks asso- 
ciated with high-mass protostars in detail. 
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Fig. 1. — Three-color images (blue, green, and red for 3.6, 4.5 and 8.0 /zm, respectively) and one dimensional 
intensity distributions (blue, green, gray and red for 3.6, 4.5, 5.8 and 8.0 /zm, respectively) for individual 
objects. In all the three colors north is up and east is right. For the these images the color contrasts were set 
with the upper and lower limits based on the specified percentage of the flux in the entire image (see text). 
The saturated regions near these positions are masked in black. The small cross in each three-color-image 
shows the position of the protostar (HPMOs, HC/UC H lis) measured using millimeter interferometry (see 
Table [2] for references) . The white curve in the three-color image indicates the region where we measured 
the flux ratios and colors. The positions where we extract the one dimensional intensity profiles are marked 
using dotted lines. In the one-dimensional profiles, the flux outside these positions is indicated using dotted 
curves. In these plots the flux is shown in logarithmic scale in MJy str -1 . 
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IRAS 05358+3543 G 192.16-3.82 




15" 30" 



Fig. 2.— Three-color images for IRAS 05358+3543, G 192.16-3.82, IRAS 20126+4104, and W 75 N with 
different contrasts from Figure [T] In all the three colors north is up and east is right. Dashed circles, curves 
and rectangles with arrow are the positions where the presence of H2 2.12 /^m emission was reported in the 
literature (Varricatt et al. 2010 for IRAS 05358+3543, G 192.16-3.82, and IRAS 20126+4104; Davis et al. 
1998 for W 75 N). The regions where we measured the flux ratios and colors are also indicated using red or 
green curves. 
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Fig. 3. — Examples of background subtraction, (left) the /4.5^m image of the objects after convolution. 
The coordinates are shown in pixels (1 pixel correspond to 0".6). The crosshair in the image shows the 
positions where we measured the one-dimensional flux distribution for background subtraction, (middle) one- 
dimensional intensity profiles at 3.6, 4.5 and 5.8 fim before subtracting the background. Gray squares show 
the range where we measured the background level and its standard deviation after background subtraction 
(i.e., the uncertainty of the flux measurement for the flux ratio maps and color-color diagrams), (right) 
same as the middle but after background subtraction. The arrows and dashed lines show the spatial ranges 
for >15er (>18.7 MJy str -1 ) for the 4.5-/xm emission. This corresponds to the region we selected for the 
color-color diagrams (see text for details). 
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Fig. 4. — The and hs/im/ h.8^,m maps for individual objects. In all the three colors north 

is up and east is right. The region is selected based on the criterion in Table [2] The range of grayscale 
for the hs/im/h^/im and h.6fj,m/l5.8fj,m maps is 0-1 for all figures but L 1527 (0-1.2 and 0-1.8 for this 
object). The cross shows the position of the protostar (HPMOs, HC/UC H lis) measured using millimeter 
interferometry (see Table [2] for references). The saturated regions near these positions are masked in white. 
Black curves are drawn in the hspm/h.&iim maps for G 35.2-0.7 N, G 192.16-3.82, IRAS 2016+4104 and 
L 1527, and these are the regions selected for color-color diagrams excluding stars near the boundary. For 
the other objects the entire region shown in the figure is used for color-color diagrams. 
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Fig. 5. — The [3.6]- [4.5] versus [3.6]- [5.8] colors measured in "Green Fuzzies" in ten high-mass star forming 
regions. The error bars are shown only for those larger than the size of the dots. The solid arrows show 
the extinction vector based on the measurements of molecular clouds ([3.6]-[4.5] and [3.6]-[5.8] of 0.09±0.03 
and 0.18±0.03, respectively, at A K > 2; Chapman et al. 2009). Gray arrows show the extinction vectors at 

±1-CT. 
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Fig. 6. — Emission deviating from the linear color correlation in the G 192.16-3.82 and IRAS 20126+4104 
regions, and their spatial distributions in the 4.5-/im image. The black and gray dots in the color-color 
diagrams show the deviating and remaining components of the [3.6]-[4.5] and [3.6]-[5.8] colors. In the 4.5- 
/im image the former is shown with thick curves, while the entire region used for the color-color diagrams 
are shown with thin curves. 
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Fig. 7. — Same as Figure [5] but with colors measured in the scattered continuum in the L 1527 outflow. 
The dashed line shows the regression line for the results of L 1527 ([3.6]-[5.8] = [3.6]-[4.5] x 1.374-0.067). 
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Observations vs. models with scattering in optically thin regions 
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Observations vs. models with optically thick, geometrically thin regions 
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Fig. 8. — Same as Figure [5] but with scattering for models of the optically thin cases (upper) and optically 
thick but geometrically thin cases (lower). The colors for the objects are shown with black dots, with error 
bars if they are larger than the dots. Gray crosses show the colors for star-disk systems (Robitallile et al. 
2006) with scattering. Open squares show the colors for a blackbody (T=4000, 1000 and 500 K from the 
lower-left to the upper-right) with scattering. The gray areas show the modeled colors, adding arbitrary 
extinction, assuming that the linear correlation observed in L 1527 is due to different extinction between 
positions. 
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Fig. 9.— Same as Figure [6] but for G 35.2-0.7 N and W 75 N regions. 
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Fig. 10. — Same as Figure [5] but with those observed in shocks associated with low-mass protostars. The 
error bars are shown only if they are larger than the dots. 



25 



H+He (W07+L99) H+He (L99) H 2 +He 




J I I I I I I I I I I I I I L 



0123 0123 0123 

[3.6]-[4.5] 

Fig. 11. — Colors of thermal H2 emission for LTE with all the transitions (solid curves) and non-LTE 
with 36/49 levels (dark gray area). The upper and lower figures are for isothermal cases and temperature 
structures determined by a power-law cooling function (A oc T s ), respectively. The left plots are for collisions 
with H+He based on the collisional rate coefficients of H and He by Wrathmall et al. (2007) and Le Boutlot 
et al. (1999), respectively; middle are the same but the coefficients of H and He both by Le Boutlot et 
al. (1999); right are the same as middle but for collisions with H2 and He. Gray solid and dashed curves 
show colors with the same density and temperature (or temperature structure determined by the power-law 
index), respectively. The light gray areas show the gap between LTE calculations for all and those with 
limited numbers of the transitions (see text). 
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Observations vs. thermal H2 (with H+He collisions) 
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Observations vs. thermal H 2 (with Hb+He collisions) 
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Fig. 12. — Same as Figure [5] but with the possible coverage for thermal H2 emission (gray area) shown in 
Figure [TT] 
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Fig. 13. — Same as Figure [5] but with diagrams for modeled thermal H2 emission for collisions with H2 and 
He shown in Figure 11 
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Fig. 14. — Same as Figure 4 but with models for PAH and fluorescent H 2 emission. See text for details of 
models. 
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Table 1 

List of Lines in the IRAC Bands 



Transition A (>m) a E„/fe (xl0 3 K ) a IRAC Band b 



— Bright H2 lines c — 



1-0 0(5) 


3.235 


7.0 


1 


2-1 0(5) 


3.438 


12.6 


1 


1-0 0(6) 


3.501 


7.6 


1 


0-0 S(14) 


3.724 


19.4 


1 


1-0 0(7) 


3.808 


8.4 


1 


0-0 S(13) 


3.846 


17.4 


1 


0-0 S(12) 


3.997 


15.5 


2 


0-0 S(ll) 


4.181 


13.7 


2 


0-0 S(10) 


4.410 


11.9 


2 


1-1 S(ll) 


4.417 


19.0 


2 


0-0 S(9) 


4.695 


10.3 


2 


1-1 S(9) 


4.953 


15.7 


2 


0-0 S(8) 


5.053 


8.7 


3 


0-0 S(7) 


5.511 


7.2 


3 


0-0 S(6) 


6.109 


5.8 


3 


0-0 S(5) 


6.909 


4.6 


4 


1-1 S(5) 


7.281 


4.6 


4 


0-0 S(4) 


8.026 


3.5 


4 


0-0 S(3) 


9.665 


2.5 


4 


— Possible atomic and ionic 


lines in 


shocks d — 


Pf 7 


3.741 


155.3 


1 


Br 7 


4.052 


151.5 


2 


Pf & 


4.654 


154.6 


2 


Fe II 


5.339 


2.7 


3 


Ni II 


6.636 


2.2 


4 


Ar II 


6.985 


2.1 


4 


Pf a 


7.460 


153.4 


4 


Ar III 


8.991 


1.6 


4 



a Based on |Draine and Bertoldi| ( |1996| l for 

H2 ,and National Institute oT Standards and 
Technology (NIST) Atomic Spectra Data Base 
{ http:/ /www. nist.gov/pml/data/asd.cfm ) for atomic 
and ionic lines. 

b 1—4 for 3 . 6, 4.5, 5.8, and 8.0 fim, respectively. Sec, e.g., 
|Reach et al.| {2006 1 | for the spectral response function. 

c The line list is based on |Smith and Rosen||2005 ) : Ybarra 
I and Lada| {2009} and calculations by jTalca'mi^tHli 
Note that the four IRAC filters cover more than 200 tran- 
sitions of molecular hydrogen, and those not listed here 
significantly contribute to the flux in the 3.6-/xm band at 
T— 3000— 4000 K, i.e. , close to the dissociation temperature 
dTakami et al.||2010| see also Section 4.2). Note that the 
1-0 S(l) line (2.12 fim) described in the text has an upper 
energy of 7.0 X 10 3 K, i.e., the same as 1-0 0(5) in the list. 

d The line list is based on [Reach et ah] |2006[ |. Note that 
the the line fluxes highly depend on the ionization state, 
which is a function of temperature, densi ty, and UV field 
produce d in high-temperature slabs (e.g., [Hollenbach and| 
|McKee|[l989t . Thus, the upper level energies do not di- 
rcctly indicate the temperature which the individual emis- 
sion lines trace. 
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Table 2: Target List 



Object 


Category 21 


R.A. (2000) b 


Dec (2000) b 


Distance 


Uncertainty of IRAC 


Region selected for color-color 












(kpc) 


fiux d 
3.6 /xm 


(MJy str" 
4.5 


5.8 /xm 


diagram 


IIRAS 05358+3543 


HPMO 


5h 39m 13.0s 


35° 


45' 51" 


1.8 


0.44 


0.53 


0.99 


>15-a" for 4.5 /xm 


IRAS 16547-4247 


HPMO 


16h 58m 17.2s 


-42° 


52' 08" 


2.9 


0.67 


0.98 


1.44 


>15-er for 4.5 /xm 


G 35.2-0.7 N 


HC HII 


18h 58m 13.0s 


1° 


40' 36" 


~2 


0.36 


0.42 


0.87 


>15-<r for 4.5 /xm 


G 192.16-3.82 


HC HII 


5h 58m 13.5s 


16° 


31' 58" 


~2 


0.14 


0.16 


0.46 


>15-(T for 4.5 /xm 


IRAS 20126+4104 


HC HII 


20h 14m 26.0s 


41° 


13' 33" 


1.7 


0.50 


0.94 


3.45 


>15-<r for 4.5 /xm, >5-<r for 5.8 /xm 


W 75 N 


UC HII 


20h 38m 36.4s 


42 o 


37' 34" 


~2 


0.28 


0.19 


0.57 


> 1.0 x 10 2 MJy str -1 for 4.5 /xm 


G 11.92-0.61 


C08 


18h 13m 58.1s 


-18° 


54' 17" 


3.8 


0.43 


1.24 


0.77 


>15-<r for 4.5 /xm 


G 298.26+0.74 


C08 


12h 11m 47.7s 


-61° 


46' 21" 




0.23 


0.22 


0.24 


>80-cr for 4.5 /xm 


G 324.72+0.34 


C08 


15h 34m 57.5s 


-55° 


27' 26" 




0.20 


0.28 


0.40 


>15-<r for 4.5 /xm 



a HPMO ... high-mass protostellar objects without any previous detection of hypercompact or ultracompact H II regions; HC/UC 
HII ... hypercompact/ulracompact H II regions; C08 ... candidates of high-mass protostars discovered in Spitzer GLIMPSE 
survey by Cyganowski et al. (2008). The categorization of HPMO, HC H II and UC H II is based on |Beuther and Shepherd| 
l(2005}. 

^Positions for HPMO, HC/UC H II are based on millimeter interferometry by|Beuther et ah] | |2002| ; |Franco-Hernandez et al.| 
l |2009^|Gibb et al.| |2003ll; |Shepherd et al.| l |2000| |200l| |2003| l; |Sollins et al.| ( |2004| |. Those for the three additional objects are 
from I Cyganows ki et al.| l |2008[ |. 

c Quoted from Cyganowski et aT| ( |2009} for G 11.92—0.61, and the above references for the others. Note that our analysis and 
discussion are independent of the distance listed here, as seen in Sections 3—6. 
d Measured at an angular resolution of ~3" (see text). 
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